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Abstract. 

We present Tf-band images of 31 radio quasars selected from the B3-VLA quasar sample. A new method has 
been developed to detect and restore the host galaxies of these quasars. We have detected the host galaxies of 16 
of them. Using these data together with previous data from the B3-VLA quasars (Carballo et al. 1998, hereafter 
Paper I) and literature data from other samples, we have analyzed a sample of 69 hosts of radio quasars, covering 
a redshift range between 0<~ z <3. 

The host galaxies are large (r c ~15kpc) and luminous elliptical galaxies (~75% of them brighter than L* K ), 
with an evolution similar to that of radio galaxies. A significant fraction (~40%) of them shows evidence of a 
possible collision/merging process. They follow a ii c —r c relation similar to that of normal elliptical galaxies. The 
morphological and photometric similarities between these galaxies and radio galaxies in this wide range of redshifts 
is a good test of the reliability of unification schemes. All of them show little evolution from 2=3 to the present 
epoch. Their K — z distribution is consistent with a no-evolution model, with a fraction of the dispersion due to 
differences in radio power. We have found a correlation between the quasar radio power and the host luminosity. 
We have also found a correlation between the host and nuclear source luminosities. These correlations could be 
induced by a physical relation between the central black hole and the bulge mass (Magorrian et al. 1998). 
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i— j , 1. Introduction Studies of samples of low luminosity and low red- 

shifts AGNs (e.g., MacKenty 1990; Gonzalez-Serrano et 

M I The studies of the host galaxies (HGs) of active galactic a l. 1993) showed a difference in the morphology of the 

■ nuclei (AGNs) can help us to answer important questions host galaxies of radio loud (elliptical) and radio quiet (spi- 

related with these peculiar objects: (i) what is the ori- ra i) sources. Preliminary studies of the HGs of quasars 

gin of the nuclear activity? (ii) what is the origin of ra- f rom groundbased telescopes apparently supported this 

dio emission? (Smith & Heckman 1990; Hutchings & Neff hypothesis to explain radio activity (e.g. Smith et al. 1986; 

1992), (iii) which AGNs can be unified by orientation ef- Veron-Cetty & woltjer 1990). These studies are limited by 

fects (standard unification schemes, Antonucci 1993; Urry the atmospheric seeing: The central point-like source con- 

& Padovani 1995) and which not? (e.g., Ellingson et al. taminates the surrounding HG, and/or even blurred by 

1991), (iv) which is the proper model to explain the it, due to the point spread function (PSF) . Due to this, 

quasars evolution? (e.g., Kauffmann & Haehnelt 2000). It these studies are very complex, especially at high z, where 

is necessary to determine the space of parameters filled by the surface brightness of the HGs drops significantly due 

the different families of AGNs in order to determine which to cosmological dimming. Observing at the near-infrared 

apparently normal galaxies could contain a dormant AGN. (NIR) reduces the central point-like contribution, since the 

ratio between the luminosity of the nuclear source and the 

~ ; ~ ~~ ~~~~ _ _ , host is minimized (Dunlop et al. 1993). This strategy has 

bend offprint requests to: b.t .Sanchez , v , , ' , 

. , „ „ . , . , . , been used by several authors (McLeod & Rieke 1994a, b; 

Visiting Astronomer, German-Spanish Astronomical J v ' ' 

Centre, Calar Alto, operated by the Max-Planck-Institute Taylor et al. 1996; Hutchings & Neff 1997; Kotilainen et al. 

for Astronomy, Heidelberg, jointly with the Spanish National 1998 ; Carballo et al. 1998; Aretxaga et al. 1998; Percival 

Commission for Astronomy. et al. 2001). At low z (z <0.4) they found that (i) the 
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Table 1. Characteristics of the data 



Date 


Telescope 


Instrument 




FOV 


Pixel Scale 


Exp. time 




seeing 


/^K lim 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


29 Sep 96 


4.2m WHT 


WHIRCAM 


~ V 


' x 1' 


0.240"/pixel 


900s 


0.03 


0.8" 


21.8 


14 Oct 97 


3.5m CAHA 


OMEGA 


~6' 


x 6' 


0.396"/pixel 


1800s 


0.06 


1.2" 


20.9 


15 Oct 97 


3.5m CAHA 


OMEGA 


~6' 


x 6' 


0.396"/pixel 


1800s 


0.16 


1.1" 


20.9 


16 Oct 97 


3.5m CAHA 


OMEGA 


~6' 


x 6' 


0.396"/pixel 


1800s 


0.12 


0.9" 


20.9 


3 Feb 99 


3.5m CAHA 


OMEGA 


-6' 


x 6' 


0.396"/pixel 


2700s 


0.09 


1.3" 


21.0 



(1) Observing date. (2) Observing telescope. (3) NIR instrument used. (4) Field of view of the images. (5) Pixel scale. (6) 
Exposure time used in the science exposures. (7) Photometric error derived from the standard deviation around the mean 
photometric zero-point per night. (8) Mean FWHM of field stars on the images along each night. (9) Mean 3a limiting surface 
brightness magnitude of the science exposures for each night. 



radio quasars inhabit large and luminous elliptical galax- 
ies (r ~10kpc and Lrgs >2L*), and (ii) a considerable 
fraction of the radio-quiet quasars inhabits also elliptical 
galaxies, and not only spiral ones, in opposition with what 
it was previously claimed. The most luminous quasars in- 
habit elliptical rather than spiral galaxies, despite their 
radio-loudness. Detailed morphological analysis at higher 
z ranges was not possible. It seems that the HGs of dif- 
ferent families of radio sources have similar morpholo- 
gies (down to z <0.4) and luminosities (down to z <1), 
which reinforces the standard unification schemes of radio 
sources. 

More detailed morphological studies have been done 
using the EST (e.g. Bahcall et al. 1994, 1995a,b, 1996, 
1997; Hutchings ct al. 1994; Disney et al. 1995; Hutchings 
& Morris 1995; McLure et al. 1999; Lehnert 1999a,b; 
Kirhakos et al. 2001; Hutchings et al. 2002). They show 
that there is a large fraction of HGs with unusual and 
irregular morphologies, undergoing a merging process or 
with evidence of such a process in their recent histories. 
However, these results were not completely conclusive due 
to: (i) the saturation and undersampling of the point-like 
source in the EST deep images (e.g., Hutchings ct al. 1995; 
McLure et al. 1999); (ii) strong contamination from emis- 
sion lines in the optical bands at certain redshifts; (iii) 
the use of bands that sample the flux down to A ^-TOOOA, 
in a range not dominated by the stellar emission, which 
limits the redshift range for studies that use optical bands 
(e.g. z <0.6 for the i?-band); (iv) studies mainly based on 
low-z samples or with a reduced number of objects. On 
the other hand, the few published studies using NICMOS 
data appear to confirm the groundbased NIR studies (e.g., 
McLeod & McLeod 2001; Kukula et al. 2001; Ridgway et 
al. 2001). 

The number of spectroscopic studies of HGs is still 
limited due mainly to the technical difficulties. The first 
attempts to determine the nature of the 'Fuzz' detected 
around the quasars were made nearly two decades ago 
(Boroson & Oke 1982, 1984, and references therein). With 
a large number of caveats, they found that the 'Fuzz' 
around QSOs has an stellar nature. More recent studies 
(Canalizo & Stockton 2000; Hughes et al. 2000; Nolan ct 
al. 2001; Canalizo & Stockton 2001; Courbin ct al. 2002) 



have confirmed the stellar nature of the extended emis- 
sion. They found that at low z these objects are dominated 
by an old stellar population (t agc >12Gyr), with possible 
traces of more recent star formation events that might be 
due to past interactions or even major mergers (Canalizo 
& Stockton 2001; Stockton & Ridgway 2001). 

There have been few attempts to study the properties 
of HGs of radio loud quasars over a large and representa- 
tive sample. Lehnert et al. (1999a) studied a sample of 43 
3CR radio quasars, using i?-band snapshots HST images. 
They detected the host for about ~ 50% of their quasars. 
These HGs present distorted morphologies, close compan- 
ions (in a ^25% of the objects), and a significant align- 
ment with radio emission. As we quoted above, i?-band 
images do not trace well the stellar population of the HGs 
at z >0.6. The use of this band combined with the depth 
of images (16 of the 19 undetected hosts are at z >1), and 
the lack of morphological classification make their results 
inconclusive. In contrast, McLure et al. (1999) and Dunlop 
et al. (2001) found, from i?-band deep HST imaging of a 
sample of z <0.4 objects, that the morphologies of the 
HGs of radio-loud and radio-quiet quasars and of radio 
galaxies were very similar (elliptical galaxies), and that 
there were no significant traces of interactions and/or dis- 
torted morphologies. Their sample covers a lower range 
of radio powers than the one studied by Lehnert ct al. 
(1999a). It might be impossible to generalize the proper- 
ties of the HGs of powerful samples of radio quasars, like 
3CR, since these objects trace only the most violent events 
of the radio source family. 

It is necessary to study the HGs of radio quasars in a 
wide range of redshifts, for complete samples with interme- 
diate radio power, trying to trace the stellar component of 
their HGs. The if-band is the best choice to compare with 
other families of radio sources since they have been previ- 
ous studied in this band (e.g., Lilly et al. 1985; Willot et 
al. 2002). We started this study in 1996 using the B3-VLA 
quasar sample (Vigotti et al. 1997; Sanchez et al. 2001). 
This sample consists of 130 radio quasars selected up to 
S408MHZ >0.1 Jy. As a part of it we published the results 
based on the study of the if-band images of 54 quasars 
of this sample (Carballo et al. 1998, hereafter Paper I), 
32 of them were analyzed in order to detect the HG and 
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determine its photometry. It was possible to recover the 
HG in 16 of them (50%); These were large and luminous 
galaxies, roughly similar to radio galaxies down to z <1. 
We present here the study of new if -band images of 31 
radio quasars from the B3-VLA sample. 

The distribution of this paper is as follows: in Sect. 2 
we present the sample and the observations and in Sect. 

3 we show the technique used to restore the HGs of the 
quasars. In Sect. 4 we present the results of this analysis. 
In Sect. 5 we discuss the evolution of HGs of radio sources, 
showing its K — z distribution, the structural parameters 
of these HGs, and its relation with radio power and nuclear 
emission. In Sect. 6 we study the relation between nuclear 
emission and radio power. We present the conclusions in 
Sect. 7. Throughout this article we have assumed a stan- 
dard cosmology with H o =50 Km s _1 Mpc -1 and g o =0.5. 
The selection of other cosmological parameters will not 
significantly affect our results 1 . 

2. Observations, data reduction and photometry 

Tabic ^ shows a summary of the observations, including 
the observing nights, the telescopes and instruments used, 
their field of view (FOV) and pixel scale, the exposure 
time, the error of the photometric calibration (<tk, ex- 
plained below), the mean seeing during each night (deter- 
mined from the mean FWHM of the stars observed along 
the night) and the mean surface brightness limit at 3er for 
each night. 

A standard NIR observing procedure was used: n un- 
registered exposures of t seconds each one were taken and 
the average image was registered. A number of k average 
images were obtained following a grid pattern around the 
object central position, with an offset of Z arcseconds. Both 
the target and the photometric standards were observed 
using the same procedure. The selection of the above- 
mentioned parameters {n,t,k and I) depended on the char- 
acteristics of the detector and on the sky background level 
of the night. For the WHIRCAM images 75 unregistered 
exposures of 2 second, and 6 average images with a 17" 
offset were taken, which allowed a total exposure time of 
900 seconds. For the October 1997 OMEGA images 25 
unregistered exposures of 2 seconds, and 9 average images 
with a 35" offset were taken. The procedure was repeated 

4 times to reach a total exposure time of 1800 s. For the 
February 1999 images the exposure time of the unregis- 
tered frames was 4 seconds and the whole procedure was 
repeated only three times. The total exposure time was 
2700 s. 

The data were reduced using standard IRAF 2 pack- 
ages. First, a dark frame, obtained before each sequence 
of exposures along the grid, was subtracted from each reg- 
istered image. For the OMEGA images, domeflats were 

1 A cosmology with H =65 Km s" 1 Mpc" 1 , fi=0.7 and 
A=0.3 will change the absolute magnitudes of the objects by 
~0.3 mag in the range of redshifts we considered 

2 IRAF is distributed by the NO AO, which is operated by 
AURA, Inc., under contract to the NSF 



obtained as described in the OMEGA observer's manual. 
For the WHIRCAM images a skyflat was built comput- 
ing the mean of all the registered images along the night, 
applying a la rejection algorithm to remove source con- 
tributions. These flat-field images were then used to cor- 
rect the pixel-to-pixel gain variations for each registered 
image. The sky background frame was then obtained for 
each sequence of k images (six for WHIRCAM, nine for 
OMEGA), computing the mean of these images, applying 
a la rejection algorithm to remove source contributions. 
The sky-flux frame was subtracted from each of the images 
of the sequence. Once sky-subtracted, the images of each 
sequence were re-centred and co-added. When more that 
one sequence was obtained for a given object, we have re- 
centred and co-added all the sequences to obtain the final 
image. 

All the images were obtained in photometric condi- 
tions. Flux calibration was carried out using UKIRT faint 
standard stars (Casali & Hawarden 1992) observed dur- 
ing each night. The number of calibration exposures was 
always more than 12 per night for the OMEGA runs, and 
4 for the WHIRCAM run (which comprised only ~2.5 
hours), taken in the course of the night at different air- 
masses. Due to the small extinction in the NIR no air mass 
correction was needed, as we did not expect significant 
differences between the derived calibration zero-points at 
different air-masses. We used the mean zero-point for each 
night in order to calibrate the exposures. The standard de- 
viation around this value yields the photometric error due 
to the calibration for each night. The if-band magnitudes 
were them obtained by standard aperture photometry, se- 
lecting for each object the aperture where the flux reaches 
the background. 

3. Analysis of the data 

As we mention in the introduction there are various prob- 
lems in the detection and restoration of the HGs of quasars 
at high z; among these we consider that the most impor- 
tant ones are: (i) the active nucleus dominates the emission 
over a wide range of wavelengths, (ii) the central source 
flux contaminates the HG flux, due to the width of the 
PSF, (iii) the angular sizes of the galaxies are smaller at 
z ~1 than at low z(<0.4), viz. in the range of a few arcsec- 
onds, and (iv) at higher redshift, the angular size increases 
again, but cosmological dimming makes the surface bright- 
ness of the HGs drop significantly following a (z + l) 4 -law 
(i.e., about 3 magnitudes from z^0 to 2~1). 

We have minimalized as much as possible the effects 
of these problems. The first problem was reduced by the 
selection of the wavelength. At A ~ 1/im the SED of the 
galaxies presents a peak, whereas the SED of the quasars 
presents a minimum (Dunlop et al. 1993). Therefore, at 
z ~1 the ratio between HG and quasar fluxes is maximal 
in the if-band (A c ff ~ 2/im). An additional advantage is 
that in the NIR the HG emission is dominated by its stellar 
component, roughly related with the mass of the galaxy. 
This is in contrast with the optical bands, where the de- 
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Fig. 1. Surface brightness profile (SBP) of the quasar B3 
0740+380C (solid circles) together with the SBP of the 
field stars scaled to the quasar intensity peak (points). 
The solid line shows the mean of the stars' SBPs, and the 
dashed line shows the mean PSF derived from fitting a 
Moffat function. 

tection of extended emission around the quasars does not 
guarantee that the HG has been detected, since there are 
other more important contributions to the extended flux, 
such as gas flow emission (e.g., Heckman et al. 1991). 

We have assumed that any extended emission detected 
around the quasars is associated with the stellar compo- 
nent of the HG. Other types of extended emission, such 
as radio-aligned components similar to those seen in radio 
galaxies (Dunlop & Peacock 1993) or synchrotron-beamed 
emission, might be detected in the if-band, although they 
are expected to contribute less than 10% to the HG flux 
(Riegler et al. 1992; Ridgway & Stockton 1997). For the 
objects at z >2 some contribution from H Q is expected, 
about ~15% (Lehnert et al. 1999b). Therefore, we assume 
that the images of the quasars at the NIR wavelengths 
consist of two components only: (i) the nuclear, point like 
source, characterized by a scaled PSF (to match its flux), 
and (ii) the host galaxy. 

Our procedure to detect and recover the HG consists 
on the following steps: 

1. We obtained for each image the surface brightness pro- 
files (SBP) of the quasar and of a number of field 
stars. We masked the possible contaminating sources 
before obtaining the SBP. This includes basically all 
the sources that are more than 3" distant from the 
quasar peak and that do not show any connection such 
as bridges or tails with the quasar. To obtain the SBP 
increases the signal to noise ratio since it implies an 
average of the brightness along the eccentric anomaly. 
We reduce the problem to one dimension, although this 
implies a loss of information. This procedure allowed 
us to study the properties of HGs that were initially 
poorly visible in the images. The SBP was obtained 
using the techniques discussed in Jedrzejewski (1987), 



implemented in the ELLIPSE task included in the sts- 
das package of IRAF. 

The PSF was modelled by a Moffat function (Moffat 
1969). The use of an analytical model allows us to give 
a better description of the PSF in the outer regions 
(normally more noisy) and let us use the PSF for a 
proper convolution with galaxy models (see below). 
To determine the best PSF model for each image we 
fitted the SBP of the field stars to the Moffat function 



'Moffat 



(r;/ ,r 1/2 ,/3) = I a 



l + (2 



V/3 



!)•( — ) 2 

r l/2 



The width of the central core is determined by r 1 / 2 , 
and [3 determines the extension of the wings. I is the 
scaling factor that determines the peak of the inten- 
sity. Once the field stars are fitted to this model we 
obtain a mean value for both r 1 ^ 2 and (3. These values 
characterize the mean PSF in the image. 
It was found that about ~10 stars per field were 
enough to determine the PSF parameters. There were 
enough stars per field in the OMEGA images, due to 
their large field of view. However, the WHIRCAM im- 
ages have only 1 or 2 stars per field; still, these were 
found to present a stable PSF image-to-image, mainly 
due to the short time period in which they were ob- 
tained (<2.5 hours). Therefore, we used all the stars 
available in all the WHIRCAM images to determine 
the PSF parameters. 

The mean PSF was used to determine the probabil- 
ity that a quasar was extended by fitting the sur- 
face brightness profile of the quasar to the function: 
PSF(r;I ) = /Moffat(r; I ,r{/ 2 , p)- This method was 
tested on 200 simulated stars, with the same charac- 
teristics of seeing, brightness and signal-to-noise as the 
observed quasars. This yields a probability distribution 
that the surface brightness of a star fits the mean PSF, 
a distribution that can be used to estimate the prob- 
ability that a quasar is extended or not. We have ap- 
plied the same test to the field stars, finding that they 
fit the mean PSF model within a 95% probability, as 
expected. We also found that when an artificial exten- 
sion is detected in a field star, it always contribute less 
than 10% to the total star flux. 

Figure ^ ihistrates this procedure. We show the SBP 
of the quasar B3 0740+380C, together with the scaled 
SBP of the field stars used to determine the mean 
PSF. There is a dispersion of the individual star SBPs 
around the mean PSF. The quasar profile is further 
from the mean PSF than the field stars, and shows a 
clear extension. This difference in the SBP compared 
with the differences found for the stars can be used to 
set a probability for the extension to be real. 
The surface brightness profile of the quasars was then 
fitted to a two-component model: a free-scale mean 
PSF model plus a galactic model. We have initially 
restricted ourselved to the two most commonly used 



Table 2. Sample of B3-VLA quasars studied: results of the analysis 



Object 


z 


S408 


14BU 

Q408 


Obs. 


K mag 


X lv 


Pext 


%ext 


Ext. 


Q 


Pc 


P d 






-^nuc 


r c 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


0006+397 


1.830 


1.15 


0.66 


Oi 


16.72+0.12 


1.98 


95.8 


19.5 


Yes 


0.55 


48 


13 


3.7 


18.5+0.2 


17.0+0.2 


2.39+0.09 


0019+431 


1.050 


2.22 


1.34 


Oi 


15.47+0.12 


0.38 


90.4 


5.5 


No 


— 


— 


— 


— 


— 


— 


— 


0022+390 


1.932 


1.10 


0.10 


Oi 


16.65+0.13 


8.70 


99.8 


41.4 


Yes 


0.38 


84 


9 


9.1 


17.6+0.2 


17.2+0.1 


2.00+0.02 


0032+423 


1.588 


0.92 


0.94 


w 


16.24+0.08 


<0.03 


90.0 


<1.0 


No 


— 


— 


— 


- 


— 


- 


— 


0110+401 


1.479 


1.08 


0.65 


Oi 


17.13+0.09 


4.48 


98.6 


51.8 


Yes 


0.40 


84 


11 


7.7 


17.9+0.2 


17.9+0.2 


3.60+0.02 


0137+401 


1.620 


0.82 


0.85 


Oi 


16.94+0.05 


0.91 


91.8 


16.5 


No 


— 


- 


- 


- 


- 


- 


- 


0143+446 


0.813 


1.71 


0.99 


Oi 


13.99+0.16 


6.89 


99.5 


11.5 


? 


0.40 


84 


11 


7.7 


16.3+0.1 


14.1+0.2 


7.09+0.02 


0144+430 


1.790 


0.92 


1.01 


w 


16.17+0.03 


<0.03 


90.0 


15.3 


No 


— 


- 


- 


- 


- 


- 


- 


0144+432 


1.260 


0.82 


0.70 


Oi 


17.12+0.18 


8.73 


99.8 


80.2 


Yes 


0.40 


84 


11 


7.7 


17.4+0.1 


18.9+0.4 


4.45+0.02 


0157+442 


0.721 


3.26 


0.91 


Oi 


16.32+0.17 


18.74 


99.9 


83.4 


Yes 


0.39 


84 


11 


7.7 


16.5+0.1 


18.3+0.4 


4.09+0.02 


0158+394 


0.780 


0.12 


0.18 


02 


16.34+0.07 


5.45 


99.1 


6.8 


No 


- 


— 


- 


- 


- 


- 


- 


0217+417 


1.430 


0.85 


0.94 


Oi 


17.18+0.15 


2.16 


96.2 


<1.0 


No 


— 


— 


— 


— 


— 


— 


— 


0219+443 


0.852 


0.89 


1.17 


w 


14.88+0.02 


2.55 


96.7 


42.3 


Yes 


0.51 


61 


13 


4.8 


15.8+0.1 


15.4+0.1 


1.73+0.15 


0226+467 


1.216 


2.45 


0.93 


Oi 


16.27+0.08 


8.23 


99.8 


44.7 


Yes 


0.37 


84 


9 


9.1 


17.2+0.2 


16.9+0.1 


5.25+0.02 


0232+411b 


0.500 


2.75 


0.99 


w 


14.27+0.02 


7.63 


99.6 


34.2 


Yes 


0.36 


84 


9 


9.1 


15.4+0.1 


14.7+0.1 


1.44+0.01 


0249+383 


1.120 


0.95 


0.34 


w 


16.17+0.06 


4.38 


98.5 


46.0 


Yes 


0.60 


40 


17 


2.3 


16.8+0.1 


17.0+0.1 


2.01+0.13 


0255+460 


1.210 


1.74 


0.77 


Oi 


17.89+0.20 


0.30 


90.2 


<1.0 


No 


— 


— 


- 


- 


- 


- 


— 


0704+384 


0.579 


2.87 


1.06 


Oi 


15.02+0.16 


17.62 


99.9 


33.5 


Yes 


0.37 


84 


9 


9.1 


16.2+0.1 


15.5+0.2 


5.91+0.02 


0729+391 


0.663 


0.26 


0.43 


Oi 


15.60+0.14 


20.38 


99.9 


91.7 


Yes 


0.57 


48 


15 


3.2 


15.7+0.1 


18.3+0.2 


3.77+0.02 


0740+380C 


1.063 


5.55 


1.29 


o 2 


15.62+0.05 


23.02 


99.9 


35.5 


Yes 


0.60 


40 


17 


2.3 


16.8+0.3 


16.1+0.2 


4.08+0.08 


0926+388 


1.630 


0.50 


1.09 


02 


17.09+0.08 


9.95 


99.9 


23.6 


Yes 


0.56 


48 


15 


3.2 


18.7+0.4 


17.4+0.1 


4.05+0.14 


1123+395 


1.470 


0.36 


0.82 


o 2 


15.72+0.05 


5.45 


99.1 


2.2 


No 


- 


- 


- 


- 


- 


- 


- 


1148+387 


1.303 


1.83 


0.94 


o 2 


14.93+0.05 


3.85 


97.9 


<1.0 


No 


— 


— 


— 


— 


— 


— 


— 


1206+439B 


1.400 


5.69 


0.85 


02 


15.10+0.05 


4.85 


98.8 


<1.0 


No 
















1315+396 


1.560 


1.16 


0.49 


02 


15.93+0.06 


4.85 


98.8 


28.3 


Yes 


0.49 


61 


13 


4.8 


17.3+0.3 


16.3+0.1 


3.20+0.11 


1435+383 


1.600 


0.54 


0.85 


02 


16.32+0.06 


4.25 


98.3 


<1.0 


No 
















2311+469 


0.745 


4.34 


0.69 


Oi 


14.79+0.06 


1.01 


92.2 


3.5 


No 
















2332+388 


3.200 


0.42 


0.94 


Oi 


17.66+0.14 


0.68 


91.0 


12.5 


No 
















2344+429 


1.556 


0.92 


0.58 


Oi 


16.30+0.07 


1.10 


92.7 


<1.0 


No 
















2349+410 


2.046 


1.38 


0.98 


Oi 


16.78+0.18 


3.18 


97.5 


13.8 


Yes 


0.48 


61 


13 


4.8 


19.1+0.4 


17.1+0.2 


3.35+0.02 


2351+456 


2.000 


2.21 


0.27 


Oi 


16.03+0.16 


6.39 


99.4 


17.9 


Yes 


0.39 


84 


11 


7.7 


18.0+0.2 


16.4+0.2 


2.54+0.08 



(1) Object name in the B3-VLA sample. (2) Object redshift (Vigotti et al. 1997). (3) Radio flux in Jy at 408 MHz. (4) Radio spectral index, aiol°, defined as S1460 = 
Swa" - "*™" ■ (5) Observing run: W, September 1996, WHT La Palma. Oi, October 1997, 3.5m Calar Alto. 2 , February 1999, 3.5m Calar Alto. (6) K-band aperture 
photometry. (7) X 2 jv of the surface brightness fitting to the PSF model. (8) percentage probability of being extended, derived from the previous \ 2 l v - (9) Contribution of 
the possible extension to the total flux. (10) Classification output (Yes=Extended, No=Not Extended). (11) Free parameter of the generalized model (a=\ for disk galaxies, 
a=0.25 for elliptical galaxies). (12) percentage probability of being an elliptical galaxy. (13) percentage probability of being a disk galaxy. (14) Ratio of both probabilities. 
(15) Host Tf-band magnitude. (16) Nucleus if-band magnitude. (17) Effective radius of the Host in arcsec. 
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analytical prescriptions for the galactic model: an ex- 
ponential Freeman (1970) law describing disk galaxies 
and a de Vaucouleurs (1948) r 1 / 4 law describing ellip- 
tical galaxies. The fitting procedure has three free pa- 
rameters: the nuclear component flux, the host galaxy 
flux, and the galaxy scale (for each analytical model). 
The galactic profile was convolved during each step of 
the fitting procedure by the mean PSF (normalized to 
one), in order to take into account the effects of the 
seeing on the shape of the profiles. In general, the nu- 
clear component has a more intense surface brightness, 
with larger signal-to-noise ratio. Therefore, we assume 
that with this method we are able to recover the flux 
of the nuclear component better than that of the ex- 
tended component. We then estimate the HG flux by 
subtracting the recovered nuclear component flux from 
the total flux. 

5. We found that the quasar profile fits both models well. 
In order to investigate this we applied the procedure to 
a sample of 1000 simulated quasars. These simulated 
quasars consist on a central point-like source, following 
a Moffat function, plus an extended source, following 
an exponential or a r 1 / 4 function. These images cover 
different ranges of the paramter space, such as different 
ratios between the central and the extended emission, 
between the effective radius and the seeing, and the 
signal to noise of the total source. They reproduce the 
basic parameters of the real data, such as background 
noise, pixel scale, and image depth. The simulations 
are explained in Sanchez (2001) and they will be pub- 
lished in a separate paper (Sanchez 2003). 

They show that it is not possible in our range of pa- 
rameters to distinguish between a spiral or an ellipti- 
cal HG by the comparison of the x 2 jv resulting from 
the profile fitting to the two above-mentioned models. 
This problem has been described before (McLure et al. 
1999). This could be due to the larger signal-to- noise 
of the central region, which dominates the fits to such 
an extent that subtle differences in the extended emis- 
sion become indistinguishable. We also found that the 
percentages of flux recovered for both the HG and the 
nuclear source are quite similar, independent of the 
model. 

6. To solve this problem we followed a different approach: 
(i) the HG surface brightness profile was obtained sub- 
tracting the mean PSF model, scaled to the nuclear 
source flux (previously obtained), and (ii) this pro- 
file was fitted to a generalized galaxy model (/ ga i oc 
e r ). This model includes the other two as particu- 
lar cases: the exponential model is recovered if a=l, 
and the r 1//4 -model if a=l/4. The values of a ob- 
tained with this method were used to classify the HGs. 
This method was applied to the simulated images, and 
yielded sharply different distributions of the a parame- 
ter for the spiral and for the elliptical galaxies (Figure 
[5J. We have used the obtained a to classify the HG 
morphologically. A similar procedure has been applied 



0.15 




0.05 



Fig. 2. Probability density functions for exponent a ob- 
tained from fits to simulated elliptical and disk galaxies. 

by McLure et al. (1999) to distingue between exponen- 
tial and r 1 / 4 models. 

The idea behind this method is that the most impor- 
tant region to distinguish between an elliptical and a 
disk galaxy is the galaxy core. The core is severely af- 
fected by the nuclear component substraction. Even in 
the case of a perfect subtraction there is an increase 
of the noise: The photon noise in this area is the sum 
of the photon noise from the galaxy and the QSO nu- 
cleus. Therefore, there is a loss of information in the 
inner part. Indeed, it is noticed in the simulations that 
the preferred a parameters of a simulated elliptical or 
spiral galaxy were not the expected values (0.25 and 
1.0, respectively) but somewhat larger ones (0.3 and 
1.15, respectively). Effects like this have been noticed 
by other authors (Meurer et al. 1995; Whitmore et al. 
1999). If we force the model to fit the data to the the- 
oretical values, both of them will fit equally badly, and 
cannot distinguish between the models. However, let- 
ting a free to be fitted we can recover the parameter 
that fits the data best. 

In Paper I we used a rather different method to detect 
and restore the HGs. Before the surface brightness profile 
of the quasars was obtained, the images were deconvolved 
by a mean PSF. The brightness profiles were then ob- 
tained and fitted to different models: a gaussian function 
(which characterizes the deconvolved point-like sources) 
and a gaussian plus galactic model (which characterizes 
a point-like plus an extended source). This method had 
to be applied due to the lower quality of the data used 
in Paper I. In order to compare the results obtanicd in 
Paper I with the present results, and to set limits to both 
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methods, we have applied both of them to the simulations 
mentioned above. 

We found that both methods are able to recover the 
flux of the HGs, although in different ranges of parame- 
ters. The method used in this article is valid for HGs for 
which the contribution to the total flux was larger than 
~9% (~12% for the method used in Paper I). In this range 
of fluxes, the magnitudes of the HGs are recovered with 
an error of about ~0.26 mag (~0.52 mag for the method 
used in Paper I). There is no tendency to over- or under- 
estimate the flux down to A'hg ~18.5 mag (~17 for Paper 
I). Moreover, with the method presented here it is possible 
to recover the effective radius of the HG, with an error of 
~25% (A parameter that we could not recover in Paper I). 
However, the current method requires (i) a good sampling 
of the PSF, using ^5-7 stars with the same brightness as 
the quasars, or a more reduced number of brighter stars, 
and (ii) that the sizes of the HGs are similar to or larger 
than the seeing, which, for quasars at z ~1 limits the use 
of this method to good-seeing images (i.e., better than 
~1.3"). The method used in Paper I, based on a decon- 
volution process, is still valid for lesser-quality images. It 
does not require a precise shampling of PSF and it can be 
used for images with seeing larger than 1.3", although it 
recovers the flux of the HG with larger uncertainties (0.52 
mag) and for brighter HGs (>12%). 

4. Results of the analysis 

Table |3 lists the sample of objects, including some basic 
parameters, such as the redshift and radio properties, and 
the A-band photometry. We have also listed the results 
from the analysis. A QSO has been classified as extended 
whenever (1) the probability that its SBP deviates from 
the mean PSF profile is larger than 95% and (2) the con- 
tribution of the extended source to the total flux is larger 
than 9%. These criteria have been selected based on the 
results from the simulations and our analysis over the field 
stars. 

We found an extension in 16 of the 31 quasars (55% of 
the sample). The extension around B3 0143+446 is consid- 
ered dubious, although it satisfies the criteria. This object 
shows a nearby star at less than 6" which could affect 
its SBP even when a proper mask has been applied. We 
will not include this object in our further analysis. Ten 
QSOs have a probability of being extended of more than 
99%. The remaining 6 quasars have probabilities between 
95% and 99%. Of theses 6, four are in the range of the 
faintest detected HGs (Ahg > 18 mag). The fraction of 
detected HGs decreases with redshift, being 75% of the 
sample at z <1 and 48% at z >1. Beyond this redshift 
the fraction of detected galaxies is roughly constant. This 
result is expected since at high redshift both the surface 
brightness and apparent luminosity of the HGs decrease, 
which makes their detection difficult. We have not set a 
magnitude limit for the undetected host galaxies, since the 
detection does not depend only on the brightness of the 
source. 



All the HGs have a larger probability of being elliptical 
galaxies than of being spirals. For an individual HG the 
probability that it is a spiral is a 7-17% and that it is 
and elliptical is a 40-84%. This is the first time that a 
morphological analysis of the HGs of a sample of radio 
quasars yields a conclusive result beyond z >0.4. Previous 
results (Kotilainen & Falomo 2000, Kotilainen et al. 1998, 
Kukula et al. 2001) directly assume a certain profile (r 1 / 4 ), 
or have inconclusive results from their analysis. Our result 
shows that the HGs of radio quasars at high redshift have 
the same morphology as the radio galaxies, supporting the 
unification schemes of both radio sources. Similar results 
were found by Taylor et al. (1996), McLure et al. (1999) 
and Dunlop et al. (2001), at low redshift (z <0.4). 

The photometry of seven of these objects was pre- 
sented in Paper I. The differences between the A-band 
magnitudes are within the errors (~0.3 mag). In four of 
them it was possible to follow a morphological analysis. 
Three have an extension detected in both studies (B3 
0704+384, B3 0740+380C and B3 1315+396). The differ- 
ences between their HG magnitudes are -0.5, -0.4 and 0.7 
(respectively), with a mean difference of 0.05+0.65 mag. 
These values are within the expected errors (~0.5 mag for 
Paper I HG magnitudes, and ~0.2 mag for current data). 
The extension of the remaining quasar (B3 0926+388) has 
been detected in the present study only, since its mag- 
nitude is much fainter than the detection limit for HGs 
reported in Paper I. 

The HG contribution to the total flux ranges between 
14% (B3 2349+410) and 92% (B3 0729+391), with a mean 
contribution of ^41%. Similar values have been reported 
in Paper I (range between 18%-83% and mean ~50%), 
and studies based on lower redshift samples for similar 
wavelengths (Taylor et al. 1996, Kotilainen et al. 1998). 
These values are lower at optical wavelengths, <20% for 
A-band images (e.g., Lehnert et al. 1999a, McLure et al. 
1999). Therefore, the extended emission is redder than the 
nuclear emission. This result is expected if a stellar compo- 
nent dominates the extended emission rather than other 
contributions (such as scattered flux from the nucleus). 
Recent results based on off-nuclear spectroscopy have con- 
firmed its stellar nature (Canalizo & Stockton 2000; Nolan 
et al. 2001; Canalizo & Stockton 2001; Courbin et al. 
2002). 

We have obtained a restored image of the HG by sub- 
tracting an image of the nuclear point-like source. This im- 
age was built assuming a Moffat function with the mean- 
PSF parameters (obtained by fitting to the field stars). Its 
flux was fixed to the value obtained from the fitting proce- 
dure for the nuclear component. The image was built us- 
ing the package ARTDATA implemented in IRAF. Figure 
13 shows the contour plots and surface brightness profiles 
for the sources with detected extension. In each panel (a) 
shows the contour plot (left) and profile (right) of the 
original source, and (b) shows the contour plot (left) and 
profile (right) of the recovered host galaxy. The solid line 
in (a) shows the surface brighteness profile of the PSF 
scaled to the quasar peak. The x 2 jv from the comparison 



8 



S.F.Sanchez and J.I. Gonzalez-Serrano: The Near-Infrared properties of the host galaxies of radio quasars 



of both profiles and the probability of being extended have 
also been included for each quasar. It is clearly seen that 
the objects present a significant extension from the com- 
parison of both profiles. The orientation of the contour 
plots is North (up) and East (left), and the field-of-view 
corresponds to 16" x 16". We have plotted the raw data be- 
fore masking for nearby companions. In some cases (e.g., 
B3 0006+397, B3 0110+401) there are close companions 
that affect the unmasked SBP. However, the extension is 
clearly seen in spatial regions not affected by these com- 
panions. This is not the case for quasar B3 0143+466, 
classified as dubious. 

Seven of the 16 objects present distortions, tails or 
bridges with possible close companions at <30" (40% of 
the sample) . These distortions could be traces of past col- 
lisions, interactions or even merging processes, which af- 
fect the host and nearby galaxies. In some of them, B3 
0219+443 and B3 0740+380C, the optical-to-NIR colours 
of the nearby objects (<30") are compatible with the 
colours of galaxies at the redshift of the quasars (Sanchez 
2001; Sanchez & Gonzalez-Serrano 2002). However, a con- 
tamination from a nearby object, either a galaxy or a faint 
star, within ~2" of the quasar could produce a similar 
effect and create an artificial distorted HG. We have es- 
timated the mean density of objects in our images using 
the number of objects detected in the February 1999 im- 
ages (the deepest images with the larger field of view). 
We found a density of 2.7+0.6 10 -3 obj/arcsec 2 , includ- 
ing both stars and galaxies. The large dispersion is due 
to the fact that these objects tend to inhabit clusters or 
groups of galaxies of different populations (see Sanchez & 
Gonzalez- Serrano 2002, also for a detailed explanation of 
the procedure to determine the density). Using these num- 
bers, we expect to find 0.043 objects within a box of 4"x4" 
around each quasars. I.e., we expect that a contaminat- 
ing source could create an artificial distorted HG in 4.3% 
of the objects. For our sample of 16 extended sources we 
do not expect any significant contribution from this effect 
(~0.64 possible contaminted objects). 

Together with data from Paper I, we have a collection 
of 29 HGs obtained from the analysis of the Jf-band im- 
ages of 60 radio quasars from the B3-VLA sample. Fifteen 
of them (Paper I) show evidence of possible interactions 
(~50% of the sample). This evidence is tenuous and should 
be confirmed by spectroscopics studies of both the host 
galaxies and the nearby companions. 

5. The properties of the HGs of radio sources 

5.1. The K — z distribution: The evolution of the HGs 

Figure 0] shows the iiT-band apparent magnitude of HGs 
of radio quasars as a function of redshift. The solid circles 
represent the HGs presented in this study, and the open 
circles the HGs presented in Paper I. The remaining points 
represent data obtained from the literature: diamonds 
from Lehnert et al. (1992), open squares from Taylor et 
al. (1996), asterisks from Kotilainen et al. (1998), stars 



from Kotilainen & Falomo (2000), triangles from Falomo 
et al. (2001) and solid squares from Kukula et al. (2001). 
The HG magnitudes from Lehnert et al. (1992) were cor- 
rected by 0.4 magnitudes to take into account the under- 
estimation of the HG flux due to the total subtraction of 
the central point-like source, following the indications of 
the authors. Magnitudes from Kotilainen et al. (1998), 
Kotilainen & Falomo (2000), Falomo et al. (2001) and 
Kukula et al. (2001) were transformed from H to K band 
using the correction suggested by authors (H — K ~0.2 
mag, Kotilainen et al. 1998). We have used only the 2 
HGs detected in the ii-band from Kukula et al. (2001) 
with a reliable result in the radial profile fitting. We have 
not considered their detections in the J-band (due to the 
uncertainties in the transformation to the i-T-band), and 
their dubious detections. This comprises a heterogeneous 
sample of 69 HGs, 29 B3-VLA HGs plus 40 HGs extracted 
from the literature. 

The shaded region shows the space occupied by the 
radio galaxies of different samples (±1.5<7 from the mean 
value), including the 3C and B2 (Class Uy) samples, stud- 
ied by Lilly et al. (1985) and Lilly (1989), the 6C and 7C 
sample, and the HzRG (High-z radio galaxies). The dis- 
tribution of the mean value is represented by the dashed 
line. All the radio galaxy data have been extracted from 
Fig. 1 of Van Breugel et al. (1999), except for the 7C 
data, obtained from Lacy et al. (2000) and Willot (2002). 
These radio galaxy samples, except for HzRG, are roughly 
complete (i.e. 3C is 98% complete). In contrast, only low-z 
quasar host samples are complete (e.g. Taylor et al. 1996). 
The quasar hosts are mainly ^50% complete, as only the 
brightest objects are detected. 

The distributions are roughly similar. Therefore, the 
HGs of different families of radio sources have similar lu- 
minosities at different cosmological epochs, down to z <3. 
However, 9 of the 70 HGs presented in this figure are 
clearly more luminous than the radio galaxies (~13% of 
the HGs sample), 6 of them at z >1. Considering also 
that only half of the HGs of the radio quasars have been 
recovered (or even less at z > 1, as we quoted above), 
we conclude that although the mean distributions seem 
to agree, the dispersion along the mean value is larger for 
the HGs. The errors of the HG magnitudes (upper-left 
box) are clearly larger than the errors of the magnitudes 
of the radio galaxies. This could account for the larger 
dispersion in the distribution in Fig. But not all the 
difference could be due to this cause: the largest errors 
were about 0.5 mag, and the difference between the mag- 
nitudes of the high-luminosity HGs and a brightest radio 
galaxies at the same redshift was ~1-1.5 mag. We need to 
explore other possibilities to explain these extreme lumi- 
nosities. 

Different evolutionary models are shown in Fig. 
0] Each solid concave line shows the evolution of 
an elliptical galaxy with different formation redshifts 
(z for =7,4.5,2.5,1.5,1.0 and 0.7, from right to left). Each 
convex line shows the evolution of a galaxy in a con- 
tinuous burst process with different formation redshifts 
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Fig. 3. Contour maps and surface brightness profiles of the quasars with detected extension. In each panel (a) is the 
contour map (left) and surface brightness (right) of the original image, and (b) is the contour map (left) and surface 
brigthness (right) of the HG image. The solid line in (a) shows the mean PSF scaled to the quasar intensity peak. 
The distances between small marks are 2" (for objects observed with OMEGA), and 1.2" (for objects observed with 
WHIRCAM). The size of the contour maps is 16"xl6" for all the images. First contour level is at 2a per pixel, with a 
separation of 0.5 mag between successive contours. 



(zfor=2.5,1.5,1.0 and 0.7, from right to left). These models 
have been included to show the most extreme situations, 
from pasive to violent evolution, and from old to young 
galaxies. The solid near-straight line shows the magni- 
tudes for a galaxy without evolution, which could be con- 
sidered the most extreme situation of a very old galaxy 
without present star formation. The evolution tracks have 
been determined using the GISSEL code (Bruzual & 
Chariot 1993), assuming an intrinsic luminosity of M* 
(Mobasher et al. 1993), and a Salpeter mass distribu- 
tion with a range of masses between 0.1 and 125 M Q 



(Salpeter 1955). The instrinsic luminosity has been se- 
lected to match the data at low-z for both families of 
HGs. We have used the E-type and Burst- type isochrone 
spectral energy distributions (ISED) from Bolzonella et 
al. (2000). 

The mean distribution of magnitudes with rcdshift is 
well described by an old E-type galaxy model, near to 
a no-evolution model. However, there is a large disper- 
sion around this mean track. The dispersion could be ex- 
plained by a range of formation redshifts, a range of differ- 
ent evolutions and/or a range of intrinsic luminosities (see 




Fig. 3. Continued 



Figure). In the first two cases, as we explained in Paper 
I, the brightest HGs would be galaxies in a violent star 
formation process. These galaxies would be blue and their 
predicted optical magnitudes would as bright as (or even 
brighter than) the observed quasar magnitudes. Therefore, 
a large range of star-formation redshifts or violent evolu- 
tion have to be ruled out as a general explanation to the 
data. In the same way, Nolan et al. (2001) and de Vrics 
et al. (2000) have presented evidence that low-z quasar 
hosts and radio galaxies are dominated by a stellar popu- 
lation of at least ~12Gyr. It appears more likely that the 
dispersion in the distribution is due to a dispersion in the 
instrinsic luminosity of the HGs. This dispersion seems 
to be larger for the HGs of radio quasars than for radio 
galaxies. 



Eales & Rawlings (1996), Rawlings (1998), and more 
recently Willot et al. (2002) found that the low radio- 
luminosity radio galaxies (6C and 7C) are fainter at the 
NIR than powerful ones (3C and B2) at the same redshift. 
This implies a correlation between the radio power and 
the NIR luminosity of the galaxy. Its effect on the K — z 
distribution is an increase in the dispersion. Therefore, 
the K — z distribution for the different families of radio 
sources is well described by a mean evolution typical for 
an old E-type galaxy (or a no-evolution model), and a 
range of absolute luminosities correlated with the radio 
power. The reason for this correlation could be a relation 
of both parameters (radio power and HG luminosity) with 
the central black hole mass (McLure et al. 1999; Dunlop 
et al. 2001). 
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Fig. 3. Continued 



This scenario describes the observed distribution accu- 
rately. Moreover, it explains the narrower distribution ob- 
served for the 3C radio galaxies due to the reduced range 
of radio luminosities (Lilly et al. 1985). In this context, we 
do not find radio galaxies as luminous in the NIR as the 
ultraluminous HGs observed here (~13% of the sample) 
due to a selection effect. However, this scenario implies a 
relation between the HG luminosities and the radio power. 

5.2. Relation between the HG luminosity and the radio 
power 

Figure [5] shows the radio power of the 29 HGs of B3-VLA 
quasars plus the 30 HGs obtained from literature with 
published flux at 1460 MHz (out of the original 40 HGs). 
The radio power has been obtained asuming a spectral 



index of a=— 1 for the radio emission, in all the objects. 
There is a clear correlation between these parameters for 
both the heterogeneous sample (r=0.81;P >99.99%), and 
the B3-VLA sample (r=0.70;P=99.88%). The B3-VLA 
subsample covers a range of radio-power/redshifts not 
very well covered before by data in the literature, extend- 
ing the studies to more powerful/higher rcdshifts objects. 

Figure shows the distribution of radio power at 
1460 MHz against the absolute magnitude of the HGs. 
These parameters are clearly correlated, both for the het- 
erogeneous sample (r=-0.74;P >99.99%) and the B3- 
VLA subsample (r=-0.56;P=99.81%). This correlation 
can be accounted for by the strong dependence on red- 
shift shown by both parameters (the absolute magni- 
tude of the host galaxies correlates with redshift, r=- 
0.60;P>99.99%). Splitting the sample by the mean red- 
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shift, wc find that there is still a correlation for the low- 
z subsample (r=-0.61; P >99.99%; n gal =37), although 
only a faint tendency for high-z subsample (r=-0.32, 
_P=83.93%; n ga i=22). It seems that there is an intrinsic re- 
lation between these two parameters, not induced by their 
depencence with redshift. Similar tendencies have been de- 
scribed for radio galaxies (Ledlow & Owen 1996, Rawlings 
et al. 1998; Lacy et al. 2000; Jarvis et al. 2001; Inskip et 
al. 2002; Willot et al. 2002), and radio-quiet quasars, of 
which the HGs at high z (z~2-3) are 2 mag fainter than 
powerful radio galaxies at the same redshift (Rawlings et 
al. 2001). 

There is a simple explanation of why we may expect 
to find that stellar luminosities of HGs are correlated with 
radio luminosity. In essence, this comes down to the fact 
that more massive objects are generally more luminous. 



In this case there are several observed correlations that 
can be used to specify the relationship between stellar 
and radio luminosity. First, the stellar luminosity of el- 
lipticals is known to correlate almost linearly with the 
central black hole mass (Magorrian et al. 1998), presum- 
ably because the stellar luminosity is very closely related 
to the stellar mass and there is a tight correlation be- 
tween host and black hole masses (Gebhart et al. 2000; 
Ferrarese & Merrit 2000). The radio luminosity is found 
to be linearly related to the narrow emission line luminos- 
ity (Willot et al. 2001) and therefore to the UV ionizing 
luminosity (Rawlings & Saunders 1991; Willot et al. 1999; 
Carballo et al. 1999). The UV ionizing luminosity is due 
to the accretion of material onto the supermassivc black 
hole. In fact, Franceschini et al. (1998) found a direct cor- 
relation between the radio luminosity and the mass of the 
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Fig. 4. 7-Gband magnitude distribution with redshift for the HGs of radio quasars: circles for the B3-VLA (solid for new data 
and open for Paper I data), diamonds for Lehnert et al. (1992), open squares for Taylor et al. (1996), asterisks for Kotilainen 
et al. (1998), stars for Kotilainen & Falomo (2000), triangles for Falomo et al. (2001) and solid squares for Kukula et al. 
(2001). The typical error of the magnitudes for the different samples is presented in the upper-left panel. The shaded region 
shows the parameter space occupied by radio galaxies of different samples (see text); the dashed line is the distribution of 
the mean values. The solid concave lines show the expected evolution of an elliptical galaxy with different formation redshifts 
(zfor=7,4.5,2.5,1.5,1.0 and 0.7, from right to left). The convex lines show the expected evolution of an galaxy with an ongoing 
burst with different formation redshifts (zf or =2.5,1.5,1.0 and 0.7, from right to left). The near straight solid line shows the 
expected magnitudes for a galaxy without evolution. 



central black hole, for a sample of low redshift galaxies. 
Therefore any correlation between stellar and radio lu- 
minosity is likely to have its cause in the fact that both 
parameters correlate positively with the black hole mass. 
McLure et al. (1999), and more recently Dunlop et al. 
(2001), reached a similar conclusion studying a sample of 
low- z HGs. Our result extends theirs on the radio-loud 
sources to a wide range of redshifts (0< z <3), increasing 
the statistical significance. 



5.3. Relation between the luminosities of the HG and 
the central source 

There are several pieces of evidence supporting a scal- 
ing of the nuclear luminosity with HG luminosity. As we 
quoted above, Magorrian et al (1998) found a correla- 
tion between the HG luminosity and the central black- 
hole mass. Since the accretion of matter is the main en- 
gine of the nuclear activity, it is expected that nuclear 
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Fig. 5. Distribution of the radio power at 1460 MHz 
against the redshift. The solid symbols are the HGs for 
the B3-VLA quasars; the open symbols are the HGs from 
the literature. 
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Fig. 6. Distribution of the radio power at 1460 MHz against 
the absolute magnitude of the HGs. The solid symbols are the 
HGs for the B3-VLA quasars; the open symbols are the HGs 
from the literature. 



luminosity scales also with the black-hole mass. Indeed, 
Dunlop et al. (2001) have shown that radio-loud quasars 
emit within a range of a 1-10% of the Eddington lumi- 
nosity. This luminosity is proportional to the black-hole 
mass. Therefore, a relation is expected between HG and 
nuclear component luminosities. On the other hand, hier- 
archical galaxy formation models of Kauffman & Hachnelt 
(2000), in which they have addressed both the formation 
of bulges/elliptical galaxies and the formation and fueling 
of their associated black holes, predict a faint correlation 
between both luminosities. 

Figure [7| shows the distribution of the absolute mag- 
nitudes of the HGs against the absolute magnitude of 
the nuclear source. There is a clear correlation for the 
heterogeneous sample (r=0.56;P>99.99%) and the B3- 



VLA subsample (r=0.57;P=99.88%). The strong depence 
that both parameters show with redshift could con- 
tribute significantly to this correlation (the absolute mag- 
nitude of the central sources also correlates with redshift, 
r=-0.89;P >99.99%). We have determined the correla- 
tion coefficient for the low and high-z subsamples, find- 
ing that there is still a correlation in both subsamples 
(r=0.45;P=99.61% and r=0.41;P=97.37%, respectively). 
In order to remove the redshift effect we have determined 
the correlation coefficients for the apparent magnitudes, 
instead of the absolute magnitudes. Figure|S]shows the dis- 
tribution of the apparent magnitudes of the HGs against 
the apparent magnitude of the nuclear source. We have 
found that the correlation is even stronger for the het- 
erogeneous sample (r=0.74;P>99.99%), although only a 
slight correlation is found for the B3-VLA subsample, 
which covers a narrower range of apparent magnitudes 
(r=0.34;P=96.79%). An intrinsic relation between both 
parameters rather than a redshift induced relation seems 
to be the correct explanation for the observed distribu- 
tions. 

The distributions in Figs [7\ and [S] could be induced by 
different biases. The QSO selection criteria as point-like 
sources could depopulate the lower-right region (bright 
HGs harboring a faint nuclear source). For the B3-VLA 
sample, the point-like criterion was applied to POSS plate 
data (Vigotti et al. 1997). Similar data were used for the 
different samples included in our heterogeneous sample 
(e.g. Drinkwater et al. 1997 for the PKS). Taking into 
account the mean B — K colours of the quasars and the 
elliptical galaxies (^2.5 and ~5 mag respectively), we do 
not expect a significant effect of this bias on the observed 
distribution. E.g., a quasar with B=17 mag, classified as 
point-like using the POSS plates could have a host galaxy 
that contributes 50% to its A"-band flux, and it would be 
undetectable on POSS plates (Phg ~20.5 mag). 

On the other hand, the upper-left region (faint HGs 
harboring a bright nuclear source) could be depopulated 
due to the incompleteness in the HG detection. Only HGs 
that contribute more than 10% to the total integrated 
flux are considered as real detections. This imposes an 
upper limit envelope on the observed distributions near 
to K ga \=K nuc +2.5. If this is the case we would expect 
an increase of the strength of the correlation with in- 
creasing incompleteness. However, as stated before, this 
strength is the same for the low-z and high-z subsamples, 
but the completeness decreases significantly (~75% and 
~40% respectively, for the B3-VLA sample). Moreover, 
this argument could only be applied to the B3-VLA data. 
The heterogeneous sample was built using HG samples 
that are mainly complete or less affected by uncomplete- 
ness (due mainly to the lower redshift range they cover). 
Indeed, the distribution of B3-VLA data in Fig. |H1 does 
not show a clear envelope at the predicted 2.5 mag from 
the K Nuc i ous =K H G line. 

These correlations confirm the tendency presented in 
previous articles for more powerful quasars to reside in 
more luminous hosts (e.g., Paper I, Kotilainen et al. 1998, 
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Fig. 7. Distribution of the K-band absolute magnitudes of the 
HGs against the absolute magnitude of the central source. The 
solid symbols are the HGs for the B3-VLA quasars; the open 
symbols are the HGs from the literature. 
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Fig. 8. Distribution of the A"-band apparent magnitudes of 
the HGs against the apparent magnitude of the central source. 
The solid symbols are the HGs for the B3-VLA quasars; the 
open symbols are the HGs from the literature. 



Kotilaincn & Falomo 2000, Taylor et al. 1996). As we 
stated above, the most probable reason for this correla- 
tion is a relation of both luminosities with central black 
hole mass. This correlation is weaker or even absent in the 
more nearby, lower luminosity and radio-quiet AGNs (e.g. 
McLeod & Rieke 1995; McLeod & Rieke 1994a,b; Percival 
et al. 2001). The smaller dynamical range of magnitudes 
covered by these samples could account for the absence of 
this correlation. 

5.4. Structural parameters of the HGs 

The mean absolute magnitude of the HGs is — -27.02±1.81 
mag, which corresponds to a luminosity of ~3L* (Mjjj-=- 
25.86 from Mobasher et al. 1993). ~75% of the ob- 



jects are brighter than L*. Their mean effective radius 
is 15.00±13.90 kpc and all but one have r >3.5 kpc. 
Capaccioli et al. (1992) found that there is a limit to the 
size of the galaxies (~4.5 kpc), below which it is rare to 
find AGNs, BCG and/or cDs. As expected, -90% of the 
HGs of our sample are above this limit. They are large 
and luminous galaxies, similar to BCGs or cDs, which are 
in the range of the brightest known galaxies. 

Elliptical galaxies show fundamental relations between 
their luminosities, their effective radii and their veloc- 
ity dispersions (Hamabe & Kormendy 1987; Mobasher 
et al. 1999). These relations define a plane in the three- 
dimensional space described above, the plane known as the 
Fundamental Plane (FP). The projection of this plane on 
the luminosity-radius plane (jj, e — r e ) determines a rela- 
tion between both parameters (Hoessel et al. 1987), with 
a slope near to ~3. This relation, of which the origin is 
still unknown, is shown by all the elliptical galaxies, and 
limits the structural models presented for these objects. 
It has been suggested that episodes of collision/merging 
could produce this relation (Capaccioli et al. 1992). 

We have determined the effective surface brightness 
of the HGs assuming a r 1 / 4 profile, correcting for cos- 
mological dimming (e.g., Kolb & Turner 1990). Figure El 
shows the distribution of the effective surface brightness 
with effective radius for the 50 HGs of our sample with 
published effective radius. There is a clear correlation be- 
tween both parameters, with a slope slightly larger than 3 
(~4.6), roughly similar to the value reported by Taylor et 
al. (1996). Our sample is biased towards bright and large 
galaxies, especially at high redshift. The lack of small and 
faint galaxies could artificially increase the slope of the 
correlation. The dashed lines show the location on this 
plane of a galaxy with constant magnitude (K=19 mag) 
and constant effective radius at different redshifts. Each 
line corresponds to a certain effective radius within the 
range 0.1"<r e <1.5". They represent the region of param- 
eters where the uncompleteness of the detection rise up to 
> 80% for the B3-VLA data (solid circles). Similar exclu- 
sion regions could be defined for each sample plotted in the 
figure. We have restricted the linear regression analysis to 
the data less affected by this bias (r c >10 kpc, n=27 HGs). 
We found a strong correlation (r=0.80, P >99.99%), of 
the form: 

fi c = 15.64 ±0 . 60 + 3.86 ±0 .46^og(r c ) 
similar to the well-known relation for the elliptical galax- 



6. Relation between the absolute magnitude of 
the nuclear source and the radio power 

Figure shows the distribution of the radio power of 
the quasars against the -ftf-band absolute magnitude of 
the nuclear source. Both parameters are clearly correlated 
for the heterogeneous sample (r=0.57;P >99.99%) and 
the B3-VLA sample (r=0.58;P =99.88%). This correla- 
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Fig. 9. Distribution of the effective surface brightness profile 
against the effective radius. The solid symbols are the HGs for 
the B3-VLA quasars; the open symbols are the HGs from the 
literature. The solid line show the best linear regression fit to 
the unbiased data. The dashed lines show the exclusion region 
for detection for the B3-VLA subsample. 



tion could be induced by the dependence of both parame- 
ters on redshift. Splitting the sample in low-z and high-z 
subsamples by the mean redshift, we found that the cor- 
relation persists for both subsamples (r=0.52,P=99.99% 
and r=0.51,P=98.20% respectively). This suggests that 
it is not induced by the redshift. The derived relation be- 
tween the if -band luminosity of the nuclear source and 
radio power is LK.qso ^u60MHz 4 - Carballo et al. (1999) 
found a similar relation between the optical-UV luminos- 
ity and radio power at 408 MHz, oc P^o^ 10 , for 
the B3-VLA quasars. These correlations imply a nearly 
one-to-one relation between the optical and NIR emission, 
which suggests that the principal mechanism of both kind 
of emissions is the same. The correlation between optical- 
NIR and radio emission in quasars has been discussed 
in Rawlings & Saunders (1991), Serjeant et al. (1998), 
Willot et al. (1999), Carballo et al. (1999) and Simpson 
& Rawlings (2000). It is outside the scope of the present 
study to discuss it in detail. 



7. Discussion and Conclusions 

We have presented in this study the fsT-band images of 
31 radio quasars from the B3-VLA sample, with a red- 
shift range between z ~0.5 and z ^3.2. A new procedure, 
based on the analysis of the surface brightness profile, has 
been proposed to estimate the flux contributions of both 
the point-like nucleus and the HG of the quasars. This 
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Fig. 10. Distribution of the radio power at 1460 MHz against 
the absolute magnitude of the nuclear source. The solid sym- 
bols are the HGs for the B3-VLA quasars; the open symbols 
are the HGs from the literature. 

procedure, tested on a wide sample of simulated images 
and field stars, makes it possible to assign a probability of 
being extended to each quasar, recover the flux and effec- 
tive radius of the HG, classify it as elliptical or spiral, and 
restore the image of the HG. The simulations and tests 
on field stars set limits to the detection and reliability of 
the recovered parameters. A different procedure, used in 
Paper I (Carballo et al. 1998), has also been tested on 
the simulations, in order to allow the use of the combined 
data. 

We have detected the HGs of 16 of the 31 quasars 
(~50% of the sample), all of them with a probability of 
being extended of higher than 95%, and 10 of more than 
99%. Together with data from Paper I, we have a sample 
of 29 HGs of B3-VLA quasars. Searching the literature, 
we have built a sample of heterogeneous HGs of radio- 
loud quasars that comprises 69 HGs. These HGs are large 
(r e ~15 kpc) and luminous (~75% of them brighter than 
L») elliptical galaxies. They follow the fj, c -r e relation for 
elliptical galaxies. About ~45% of the B3-VLA HGs show 
distortions and tails in their images that could be evidence 
of a merging/colliding process. This result agrees with the 
hypothesis that considers a merging process as the origin 
or triggering mechanism of nuclear activity (e.g. Smith 
& Heckman 1990; Hutching & Neff 1992). There is other 
evidence that supports this scenario: (i) the quasars are 
located in overdensities of galaxies compatible with being 
at the same redshifts (e.g. Sanchez & Gonzalez-Serrano 
1999; Sanchez 2001), and in a spatial position where a 
collision is more likely to produce a merger (Sanchez & 
Gonzalez-Serrano 2002); (ii) the traces of recent star for- 
mation processes induced by interactions (e.g. Nolan et al 
2001, Canalizo & Stockton 2000); (iii) the need of a vio- 
lent process that produces an infall of mass in the inner 
region to explain nuclear activity (e.g. Bekki & Noguchi 
1994). 
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However, it is not clear that the fraction found could 
be considered as a significant excess of merger galax- 
ies, even if all of them were confirmed. The fraction of 
merger/irregular galaxies in the local universe is ~7%, 
based on classical studies like the Shapley-Ames Catalog 
(Shapley & Ames 1932). But this fraction is known to in- 
crease with redshift. E.g, van den Bergh et al. (1996) and 
Lee & Hwang (2000) found that ^39% of the galaxies on 
the HDF-N and HDF-S show irregular/merger morpholo- 
gies. We do not have a proper sample of inactive galaxies 
with a similar redshift range and detailed morphological 
classification to known if the fraction of mergers found is 
or IS not a significant excess. 

All the HGs are elliptical galaxies morphologically sim- 
ilar to radio galaxies. The mean luminosity evolution along 
a wide range of rcdshifts (0< z <3) is also similar to the 
evolution shown by radio galaxies. This evolution is well 
described by a model with a single burst of star formation 
at very high redshift (zf OT >7) and little or no evolution 
since that period. It seems that the galaxies have been 
formed beyond z ~3, without major burst process down 
to this redshift. This supports the unification schemes of 
radio sources, since both families of radio sources are sim- 
ilar in features that do not depend on the angle-of-view. 

The HGs present a larger dispersion in mean evolution 
compared with the dispersion found in radio galaxies, es- 
pecially beyond z >1. This dispersion is most probably 
due to a correlation between the absolute magnitude of 
the HGs and the radio power. Therefore, the most pow- 
erful radio quasars inhabit the most luminous HGs. The 
absolute magnitudes of the host and nuclear source (and 
this last parameter and radio power) are also correlated. 
These correlations are most probably due to a correlation 
of all three parameters with the central black hole mass 
(McLure et al 1999). Dunlop et al. (2001) found that this 
tendency is seen both in radio-loud and radio-quiet sam- 
ples, which could imply that the same mechanism pro- 
duces the radio emission in both kinds of AGN. 

In a toy scenario, succesive merging processes could 
increase the mass of the spheroidal component of the HG, 
increasing its luminosity. These merging processes would 
feed the nuclear AGN, increasing also the mass of the cen- 
tral black hole. This will also increase the nuclear lumi- 
nosity and radio power. This scenario is consistent with 
the result that indicates that radio-quiet quasars inhabit 
less-populated areas than radio-loud quasars (Hall et al. 
1998). Sanchez & Gonzalez-Serrano (1999) and Sanchez 
& Gonzalez- Serrano (2002) have found that QSOs inhabit 
the outer regions of clusters (instead of the core) where a 
collision is more likely to produce a merger. Recent mor- 
phological analysis of HGs of different families of AGNs at 
z<0.3 found no traces of recent merging processes (Dunlop 
et al. 2001). This could indicate an evolution in the merg- 
ing rate, although much more detailed studies are needed 
before a conclusion can be reached. 

The global scenario suggested by the results presented 
in this article has to be tested. The first obvious step would 
be to increase the fraction of detected HGs in order to re- 



duce the biases, and to extend the studies to other samples 
with higher and lower radio powers. Photometric studies 
of the HGs at other wavelengths, and spectroscopy of them 
and their near companions are needed to determine bet- 
ter the evolution of these galaxies, and search for possible 
traces of merging processes. 
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